Ground penetrating radar and electrical resistivity tomography surveys with a subsequent intrusive investigation in search for the missing Beaumont children in Adelaide, South Australia.

Jane, Arnna, and Grant Beaumont went missing from Glenelg Beach in Adelaide, South Australia on 26 January (Australia Day) 1966. Despite multiple land and sea searches over nearly 60 years, the children have not been found. New credible eyewitness testimony led to a site of interest at the now disused New Castalloy factory in North Plympton, Adelaide. This site has a complex stratigraphy of anthropogenic fill, which made ground penetrating radar (GPR) investigations unpromising. Electrical resistivity tomography (ERT), while not commonly used in a forensic capacity, provided an alternative approach that allowed suitable depth penetration to resolve a feature of interest, which was subsequently excavated by the South Australia Police. This feature did contain organic, and animal remains but, sadly, not the grave of Jane, Arnna, and Grant Beaumont. However, this investigation highlights the potential to use ERT in a forensic capacity, as well as the limitations of using geophysical techniques for covert burial detection.

Monitoring of simulated clandestine graves of victims using UAVs, GPR, electrical tomography and conductivity over 4-8 years post-burial to aid forensic search investigators in Colombia, South America.

In Colombia there are estimated to be over 121,000 missing people and victims of forced disappearances. Forensic investigators therefore need assistance in determining optimal detection techniques for buried victims, to give victims' families closure and for the wider community to see that justice is being served. Previous research has created 12 controlled simulated clandestine graves of typical Colombian murder victim scenarios at 0.5 m - 1.2 m depths in savannah and rainforest sites in Colombia. The 0-3 years of geophysical monitoring results of were published, with this paper reporting on 4-8 years monitoring of both UAV drone results and geophysical data. The UAV results from the year 8 survey, published for first time from Colombia, showed that the simulated graves could still be located using NDVI and NIR multi-spectral data, but not using optical or other multi-spectral data. The 0-3 years of geophysical data found the simulated clandestine graves could be detected with electrical resistivity and GPR methods, with the 4-8 year surveys evidencing that they could still be detected using bulk ground conductivity surveys, GPR horizontal time slice datasets and 2D ERT profiles. Research implications suggest initial use of UAV remote sensing technology to pinpoint likely search areas, before subsequent ground reconnaissance, geophysical surveys and their interpretation, before intrusive investigation methods are employed for detecting missing and disappeared persons in Colombia.

The suitability of using domestic pigs (Sus spp.) as human proxies in the geophysical detection of clandestine graves.

Research in many forensic science fields commonly uses domestic pigs (Sus spp.) as proxies for human remains, due to their physiological and anatomical similarities, as well as being more readily available. Unfortunately, previous research, especially that which compares the decompositional process, has shown that pigs are not appropriate proxies for humans. To date, there has not been any published research that specifically addresses whether domestic pigs are adequate human proxies for the geophysical detection of clandestine graves. As such, the aim of this paper was to compare the geophysical responses of pig cadavers and human donor graves, in order to determine if pigs can indeed be used as adequate human proxies. To accomplish this, ground penetrating radar (GPR) and electrical resistivity tomography (ERT) responses on single and multiple pig cadaver graves were compared to single and multiple human donor graves, all of which are in known locations within the same geological environment. The results showed that under field conditions, both GPR and ERT were successful at observing human and pig burials, with no obvious differences between the detected geophysical responses. The results also showed that there were no differences in the geophysical responses of those who were clothed and unclothed. The similarity of the responses may reflect that the geophysical techniques can detect graves despite what their contents are. The study implications suggest that experimental studies in other soil and climate conditions can be easily replicated, benefiting law enforcement with missing persons cases.

Influence assessment of urban expansion on groundwater level fluctuations in Gandhinagar, Gujarat, India.

Being the state capital of Gujarat, Gandhinagar is snowballing urban population, resulting in overexploitation of groundwater resources and consequent decline in local groundwater level. The key objective of the current research is to understand the impact of urban expansion on the groundwater level of Gandhinagar district for the last 3 decades. Long-term land use/land cover (LULC) alterations using Landsat images (1991-2021) reveal a 234% increase in overall built-up area and it is more prominent in western and southern parts than the eastern part of study area till 2021 due to urban sprawl of adjacent Ahmedabad City. Spatial distribution of groundwater levels exhibits the same pattern of groundwater level drop as that of the urban expansion and the drop of maximum depth of the groundwater level has also observed during study tenure. Rapid population growth indicates inevitable urban densification which may lead to increase in groundwater abstraction and consequent groundwater level depletion of Gandhinagar City in near future. The scenario may be worsened due to the reduction in groundwater recharge area owing to enhancement of impervious surfaces. A negative correlation is established between groundwater level and respective built-up areas from 1991 to 2011. After 2001, groundwater levels in some areas showed a rising trend and the number of those locations have increased from 2001 to 2021, indicating a sufficient supply of surface water, meeting the escalating water demand and subsequent reduction in groundwater abstraction. High fluoride content was found in many groundwater samples collected from Gandhinagar's shallow unconfined aquifer. In lieu of almost unperturbed natural groundwater recharge, built-up expansion, rising population, and over-exploitation result in groundwater level depletion in both shallow and deeper aquifers. To replenish the already depleted groundwater level and for sustainable water supply, an integration of rainwater-surface water-groundwater management plan and sustainable urban management plan is highly required. The future sustainable urban-groundwater management plan of Gandhinagar City must emphasized on the expansion of green and permeable spaces for groundwater recharge, mandatory rainwater harvesting system in every building possible, suitable area demarcation for artificial recharge, and identification of areas which are less prone to groundwater level depletion for city expansion. The outcomes of the present study will help the decision-makers to prepare inclusive and resilient urban management plan to accomplish the 6th and 11th Sustainable Development Goals of United Nations by 2030.

An integrated groundwater resource management approach for sustainable development in a tropical river basin, southern India.

Evaluation of aquifer potential is essential, as the potable water demand has increased globally over the last few decades. The present study delineated different zones of groundwater potential and groundwater quality of the Kallada River basin (KRB) in southern India, using geo-environmental and hydrogeochemical parameters, respectively. Geo-environmental variables considered include relative relief, land use/land cover, drainage density, slope angle, geomorphology, and geology, while hydrogeochemical parameters include pH, electrical conductivity (EC), Cl-, Fe3+, and Al3+ concentrations. Analytical hierarchy process (AHP) was used for categorizing groundwater potential and quality zones. Nearly 50% of KRB is categorized as very high and high groundwater potential zones, occupying the western and midland regions. The central and west-central parts of KRB are characterized by excellent groundwater quality zones, while the eastern and western parts are characterized by good and poor groundwater quality zones, respectively. By integrating the groundwater potential and groundwater quality, sustainable groundwater management is observed to be necessary at about 54% of the basin, where site-specific groundwater management structures such as percolation ponds, injection wells, and roof water harvesting have been proposed using a rule-based approach. This integrated groundwater potential-groundwater quality approach helps policymakers to implement the most suitable management strategies with maximum performance.

Mapping, intensities and future prediction of land use/land cover dynamics using google earth engine and CA- artificial neural network model.

Mapping of land use/ land cover (LULC) dynamics has gained significant attention in the past decades. This is due to the role played by LULC change in assessing climate, various ecosystem functions, natural resource activities and livelihoods in general. In Gedaref landscape of Eastern Sudan, there is limited or no knowledge of LULC structure and size, degree of change, transition, intensity and future outlook. Therefore, the aims of the current study were to (1) evaluate LULC changes in the Gedaref state, Sudan for the past thirty years (1988-2018) using Landsat imageries and the random forest classifier, (2) determine the underlying dynamics that caused the changes in the landscape structure using intensity analysis, and (3) predict future LULC outlook for the years 2028 and 2048 using cellular automata-artificial neural network (CA-ANN). The results exhibited drastic LULC dynamics driven mainly by cropland and settlement expansions, which increased by 13.92% and 319.61%, respectively, between 1988 and 2018. In contrast, forest and grassland declined by 56.47% and 56.23%, respectively. Moreover, the study shows that the gains in cropland coverage in Gedaref state over the studied period were at the expense of grassland and forest acreage, whereas the gains in settlements partially targeted cropland. Future LULC predictions showed a slight increase in cropland area from 89.59% to 90.43% and a considerable decrease in forest area (0.47% to 0.41%) between 2018 and 2048. Our findings provide reliable information on LULC patterns in Gedaref region that could be used for designing land use and environmental conservation frameworks for monitoring crop produce and grassland condition. In addition, the result could help in managing other natural resources and mitigating landscape fragmentation and degradation.

Constraining the Rosalind Franklin Rover/Ma_MISS Instrument Capability in the Detection of Organics.

The Mars Multispectral Imager for Subsurface Studies (Ma_MISS) instrument is a miniaturized visible and near-infrared spectrometer that is integrated into the drilling system of the ESA Rosalind Franklin rover, which is devoted to subsurface exploration on Mars. Ma_MISS will acquire spectral data on the Martian subsurface from excavated borehole walls. The spectral data collected by Ma_MISS on unexposed rocks will be crucial for determination of the composition of subsurface rocks and optical and physical properties of materials (i.e., grain size). Ma_MISS will further contribute to a reconstruction of the stratigraphic column and acquire data on subsurface geological processes. Ma_MISS data may also inform with regard to the presence of potential biomarkers in the subsurface, given the presence of organic matter that may affect some spectral parameters. In this framework, we performed a wide range of measurements using the laboratory model of the Ma_MISS to investigate mineral/organic mixtures in different proportions. We prepared mixtures by combining kaolinite and nontronite with glycine, asphaltite, polyoxymethylene, and benzoic acid. These organic compounds show different spectral characteristics in the visible and near-infrared; therefore their presence can be detected by the Ma_MISS instrument. Our results indicate that the Ma_MISS instrument can detect organic material down to abundances of around 1 wt %. In particular, the data collected on low-concentration mixtures show that, by analyzing sediments with a grain size smaller than the Ma_MISS spatial resolution, the instrument can still discern those points where organic matter is present from points with exclusive mineral composition. The results also show that a collection of multiple contiguous measurements on a hypothetical borehole wall could help indicate the presence of organic matter in clay-rich soils if present.

Linking health to geology-a new assessment and zoning model based on the frame of medical geology.

With the growing concerns about the Earth's environment and human health, there has been a surge in research focused on the intersection of health and geology. This study quantitatively assesses the relationship between human health and geological factors using a new framework. The framework considers four key geological environment indicators related to health: soil, water, geological landform, and atmosphere. Results indicate that the atmospheric and water resource indicators in the study area were generally favorable, while the scores of geological landforms varied based on topography. The study also found that the selenium content in the soil greatly exceeded the local background value. Our research underscores the importance of geological factors on human health, establishes a new health-geological assessment model, and provides a scientific foundation for local spatial planning, water resource development, and land resource management. However, due to varying geological conditions worldwide, the framework and indicators for health geology may need to be adjusted accordingly.

A new biologic paleoaltimetry indicating Late Miocene rapid uplift of northern Tibet Plateau.

The uplift of the Tibet Plateau (TP) during the Miocene is crucial to understanding the evolution of Asian monsoon regimes and alpine biodiversity. However, the northern Tibet Plateau (NTP) remains poorly investigated. We use pollen records of montane conifers (Tsuga, Podocarpus, Abies, and Picea) as a new paleoaltimetry to construct two parallel midrange paleoelevation sequences in the NTP at 1332 ± 189 m and 433 ± 189 m, respectively, during the Middle Miocene [~15 million years ago (Ma)]. Both midranges increased rapidly to 3685 ± 87 m in the Late Miocene (~11 Ma) in the east, and to 3589 ± 62 m at ~7 Ma in the west. Our estimated rises in the east and west parts of the NTP during 15 to 7 Ma, together with data from other TP regions, indicate that during the Late Miocene the entire plateau may have reached a high elevation close to that of today, with consequent impacts on atmospheric precipitation and alpine biodiversity.

Radical change in tectonic regime and paleo-environment at the end of Neoarchean evidenced by a unique metal co-existing deposit.

Banded iron formation and Cu-Zn sulfide deposits within volcanic-argillaceous sequences (as volcanogenic massive sulphide deposits (VMS)-like type) occur together in the Qingyuan greenstone belt of the North China Craton, recording the first appearance of oxidized ores and sulfide ores co-existing in the early Earth. The unique metal co-existing deposits should meet two requirements: tectonic setting and sedimentary environment. As regards to tectonic setting, plate-like tectonics might have started since the end of the Neoarchean because continents had grown large enough and there occurred volcanic arcs and backarc basins similar to modern ones in a way. Partial melting of subducted continental crust is conductive to providing ore-forming elements. As for sedimentary environment, late Neoarchean seawater was rich in Fe2+ and anoxic. Instantaneous oxidation of the seawater resulted possibly from frequent submarine volcanic eruptions and facilitated precipitation of the banded iron formation. At this point, it is also favorable for the enrichment of Cu and Zn ions in seawater. The VMS-like deposits tended to form when the seawater was reduced again. Studies of isotopic elements like sulfur, oxygen, iron and silicon support the above geological processes. It is shown that the geologic conditions only existed in the late Neoarchean and Paleoproterozoic for a short period of time. The banded iron formations disappeared around 1.85 Ga, and the associated sulfide metal deposits also became dominant sedimentary exhalative deposits in the meso-Neoproterozoic Boring Billion, as a result of increasing oxidation of the oceans and the increasing maturity of the continental crust. This study is significant not only for decoding metallogenic genesis but also helping understand rapid change in Precambrian tectonic regimes and Earth's environments.

Seismic multi-hazard and impact estimation via causal inference from satellite imagery.

Rapid post-earthquake reconnaissance is important for emergency responses and rehabilitation by providing accurate and timely information about secondary hazards and impacts, including landslide, liquefaction, and building damage. Despite the extensive collection of geospatial data and satellite images, existing physics-based and data-driven methods suffer from low estimation performance due to the complex and event-specific causal dependencies underlying the cascading processes of earthquake-triggered hazards and impacts. Herein, we present a rapid seismic multi-hazard and impact estimation system that leverages advanced statistical causal inference and remote sensing techniques. The unique feature of this system is that it provides accurate and high-resolution estimations on a regional scale by jointly inferring multiple hazards and building damage from satellite images through modeling their causal dependencies. We evaluate our system on multiple seismic events from diverse countries around the globe. Our results corroborate that incorporating causal dependencies significantly improves large-scale estimation accuracy for multiple hazards and impacts compared to existing systems. The results also reveal quantitative causal mechanisms among earthquake-triggered multi-hazard and impact for multiple seismic events. Our system establishes a new way to extract and utilize the complex interactions of multiple hazards and impacts for effective disaster responses and advancing understanding of seismic geological processes.

Plate motion and a dipolar geomagnetic field at 3.25 Ga.

The paleomagnetic record is an archive of Earth's geophysical history, informing reconstructions of ancient plate motions and probing the core via the geodynamo. We report a robust 3.25-billion-year-old (Ga) paleomagnetic pole from the East Pilbara Craton, Western Australia. Together with previous results from the East Pilbara between 3.34 and 3.18 Ga, this pole enables the oldest reconstruction of time-resolved lithospheric motions, documenting 160 My of both latitudinal drift and rotation at rates of at least 0.55°/My. Motions of this style, rate, and duration are difficult to reconcile with true polar wander or stagnant-lid geodynamics, arguing strongly for mobile-lid geodynamics by 3.25 Ga. Additionally, this pole includes the oldest documented geomagnetic reversal, reflecting a stably dipolar, core-generated Archean dynamo.

Mapping global kimberlite potential from reconstructions of mantle flow over the past billion years.

Kimberlites are the primary source of economic grade diamonds. Their geologically rapid eruptions preferentially occur near or through thick and ancient continental lithosphere. Studies combining tomographic models with tectonic reconstructions and kimberlite emplacement ages and locations have revealed spatial correlations between large low shear velocity provinces in the lowermost mantle and reconstructed global kimberlite eruption locations over the last 320 Myr. These spatial correlations assume that the lowermost mantle structure has not changed over time, which is at odds with mantle flow models that show basal thermochemical structures to be mobile features shaped by cold sinking oceanic lithosphere. Here we investigate the match to the global kimberlite record of stationary seismically slow basal mantle structures (as imaged through tomographic modelling) and mobile hot basal structures (as predicted by reconstructions of mantle flow over the past billion years). We refer to these structures as "basal mantle structures" and consider their intersection with reconstructed thick or ancient lithosphere to represent areas with a high potential for past eruptions of kimberlites, and therefore areas of potential interest for diamond exploration. We use the distance between reconstructed kimberlite eruption locations and kimberlite potential maps as an indicator of model success, and we find that mobile lowermost mantle structures are as close to reconstructed kimberlites as stationary ones. Additionally, we find that mobile lowermost mantle structures better fit major kimberlitic events, such as the South African kimberlite bloom around 100 Ma. Mobile basal structures are therefore consistent with both solid Earth dynamics and with the kimberlite record.

Geophysical imaging of buried human remains in simulated mass and single graves: Experiment design and results from pre-burial to six months after burial.

In this study, we present an experiment design and assess the capability of multiple geophysical techniques to image buried human remains in mass and individual graves using human cadavers willingly donated for scientific research. The study is part of a novel, interdisciplinary mass grave experiment established in May 2021 which consists of a mass grave with 6 human remains, 3 individual graves and 2 empty control graves dug to the same size as the mass grave and individual graves. Prior to establishing the graves, we conducted background measurements of electrical resistivity tomography (ERT), electromagnetics (EM), and ground penetrating radar (GPR) while soil profiles were analyzed in situ after excavating the graves. All graves were also instrumented with soil sensors for monitoring temporal changes in soil moisture, temperature, and electrical conductivity in situ. Measurements of ERT, EM and GPR were repeated 3, 37, 71 and 185 days after burial with further repeated measurements planned for another twelve months. ERT results show an initial increase in resistivity in all graves including the control graves at 3 days after burial and a continuous decrease thereafter in the mass and individual graves with the strongest decrease in the mass grave. Conductivity distribution from the EM shows a similar trend to the ERT with an initial decrease in the first 3 days after burial. Distortion in linear reflectors, presence of small hyperbolas and isolated strong amplitude reflectors in the GPR profiles across the graves is associated with known locations of the graves. These initial results validate the capability of geoelectrical methods in detecting anomalies associated with disturbed ground and human decay while GPR though show some success is limited by the geology of the site.

Giant paleo-seafloor craters and mass wasting associated with magma-induced uplift of the upper crust.

Giant seafloor craters are known along many a continental margin with recurrent mass-wasting deposits. However, the impact of breakup-related magmatism on the evolution of such craters is barely understood. Using high-quality geophysical datasets, this work examines the genetic relationship among the location of magmatic sills, forced folds and the formation of giant paleo-seafloor craters underneath an ancient mass-transport complex in the Møre and Vøring basins, offshore Norway. The data reveal that forced folding of near-seafloor strata occurred because of the intrusion of several interconnected magmatic sills. Estimates of 1-dimensional uplift based on well data show that uplift occurred due to the intrusion of magma in Upper Cretaceous to Lower Eocene strata. Our findings also prove that subsurface fluid plumbing associated with the magmatic sills was prolonged in time and led to the development of several vertical fluid flow conduits, some of which triggered mass wasting in Neogene to Recent times. The repeated vertical expulsion of subsurface fluids weakened the strata on the continental slope, thereby promoting mass wasting, the selective cannibalization of the paleo-seafloor, and the formation of elongated craters at the basal shear zone of the mass-transport complex. Significantly, the model presented here proves a close link between subsurface magmatic plumbing systems and mass wasting on continental margins.

Aeolian sediment transport rates in the middle reaches of the Yarlung Zangbo River, Tibet Plateau.

Aeolian sediment emission from surfaces and subsequent transport are important geological processes. The Tibet Plateau experiences strong aeolian activity in areas such as the Yarlung Zangbo River basin. The dust storms have caused grounding of aircraft, highway closures, and other consequences for the region's residents. However, few researchers have studied this activity, which means that little knowledge is available on aeolian activity to support efforts to mitigate or prevent aeolian disasters. We measured aeolian sediment transport in the middle reaches of the Yarlung Zangbo River from 2020 to 2021. Field observations showed spatial and temporal variation of the sediment transport rate, with the greatest aeolian sediment transport in spring and winter. The largest total aeolian sediment transport rate occurred over sandy desert, with the smallest emission by a floodplain grassland. The change in sediment transport rate with height followed an exponential function, but the coefficients differed among landscapes. The mean sediment transport rate was greatest above shifting sand near riverbanks (0.21 kg m-1 d-1), where the sand is exposed in the winter and spring, followed by shifting floodplain sands (0.13 kg m-1 d-1), and was lowest above a floodplain grassland (0.03 kg m-1 d-1). Mean grain size also decreased with increasing height above 0.25 m, with a minimum mean grain size (about 52.6 µm) at 3.0 m above a floodplain grassland, and maximum mean grain size (about 100.2 µm) at 3.0 m above a floodplain shifting sand surface. The spatial variation in sediment transport rates and grain size related to the proportion of fine particles in the surface material. By comparing the aeolian sediment transport over different landscapes, we found that river banks and floodplains, which had rich deposits of very fine sand, silt, and clay, were the major sources of dust in this region. Our results indicate that efforts to mitigate or prevent aeolian disasters require a focus on riverbank and floodplain deposits.

Fine Complex Geological Structure Interpretation Based on Multiscale Seismic Dip Constraint.

With the development of seismic exploration technology, geological structure interpretation has become more and more refined, whereas random noise interference, subsalt weak seismic reflection signals, and other issues are also gradually emerged at the same time, which resulted in traditional geological structure interpretation accuracy reduction only relying on single seismic data. A novel technical process integrated time-frequency decomposition of seismic data, seismic dip constraints, and geological structure interpretation is proposed in this paper which is named multiscale seismic dip constraint geological structure interpretation. The technical process contains five steps which first use the basis tracking spectrum decomposition technology to convert the seismic data into the time-frequency domain and then decompose the raw seismic data into coarse scale, fine scale, and deliberate scale through window and threshold methods. Subsequently, execute local layer dip calculation with Hilbert transform and geological structure interpretation on seismic data of different scale, respectively. At last, perform geological structure attribute fusion to obtain fine geological structure interpretation. Synthetic data test and field data test show that through multiscale time-frequency decomposition, high-frequency noise interference can be removed and the subsalt seismic weak signal can be enhanced, and then, high-precision fine complex geological structure interpretation can be obtained with seismic dip constraint. Therefore, the technical process proposed in this paper is effective and can be widely applied in the interpretation of field seismic data.

Monitoring of simulated clandestine graves of dismembered victims using UAVs, electrical tomography, and GPR over one year to aid investigations of human rights violations in Colombia, South America.

In most Latin American countries, there are significant numbers of missing people and forced disappearances, over 120,000 in Colombia alone. Successful detection of shallow buried human remains by forensic search teams is difficult in varying terrain and climates. Previous research has created controlled simulated clandestine graves of murder victims to optimize search techniques and methodologies. This paper reports on a study on controlled test site results over four simulated dismembered victims' clandestine graves as this is sadly a common scenario encountered in Latin America. Multispectral images were collected once post-burial, electrical resistivity surveys were collected 4 times, and ground penetrating radar (GPR) surveys were collected three times up to the end of the 371 day survey monitoring period. After data processing, results showed that the multispectral data set could detect the simulated clandestine and control graves, with electrical resistivity imaging relative high resistances over some of the simulated graves but not over the empty control graves. GPR results showed good imaging on the Day 8 surveys, medium imaging on the Day 294 surveys, and medium to good imaging on the Day 371 surveys. Study implications suggest that, while clandestine graves of dismembered homicide victims would likely result in smaller-sized graves when compared to graves containing intact bodies, these graves can still potentially be detected using remote sensing and geophysical methods.

Numerical modeling as a tool for evaluating the renewability of geothermal resources: the case study of the Euganean Geothermal System (NE Italy).

Renewable natural resources are strategic for reducing greenhouse gas emissions and the human footprint. The renewability of these resources is a crucial aspect that should be evaluated in utilization of scenario planning. The renewability of geothermal resources is strictly related to the physical and geological processes that favor water circulation and heating. In the Veneto region (NE Italy), thermal waters of the Euganean Geothermal System are the most profitable regional geothermal resource, and its renewability assessment entails the evaluation of fluid and heat recharge, regional and local geological settings, and physical processes controlling system development. This renewability assessment is aimed at defining both the importance of such components and the resource amount that can be exploited without compromising its future preservation. In the second part of the twentieth century, the Euganean thermal resource was threatened by severe overexploitation that caused a sharp decrease in the potentiometric level of the thermal aquifers. Consequently, regulation for their exploitation is required. In this work, the renewability of the Euganean Geothermal System was assessed using the results from numerical simulations of fluid flow and heat transport. The simulations were based on a detailed hydrogeological reconstruction that reproduced major regional geological heterogeneities through a 3D unstructured mesh, while a heterogeneous permeability field was used to reproduce the local fracturing of the thermal aquifers. The model results highlight the role played by the resolved structural elements, in particular the subsurface high-angle faults of the exploitation field, and by the anomalous regional crustal heat flow affecting the central Veneto region.